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Abstract. We discuss the results of two XMM-Newton observations of the open cluster NGC 6604 obtained in April and 
September 2002. We concentrate mainly on the multiple system HD 167971 (05-8V + 05-8V + (081)). The soft part of 
the EPIC spectrum of this system is thermal with typical temperatures of about 2 x 10® to 9 x 10® K. The nature (thermal 
vs non-thermal) of the hard part of the spectrum is not unambiguously revealed by our data. If the emission is thermal, the 
high temperature of the plasma (^ 2.3 x 10^ to 4.6 x 10^ K) would be typical of what should be expected from a wind-wind 
interaction zone within a long period binary system. This emission could arise from an interaction between the combined winds 
of the 05-8V + 05-8V close binary system and that of the more distant 081 companion. Assuming instead that the hard part of 
the spectrum is non-thermal, the photon index would be rather steep (~ 3). Moreover, a marginal variability between our two 
XMM-Newton pointings could be attributed to an eclipse of the 05-8V -l- 05-8V system. The overall X-ray luminosity points 
to a significant X-ray luminosity excess of about a factor 4 possibly due to colliding winds. Considering HD 167971 along with 
several recent X-ray and radio observations, we propose that the simultaneous observation of non-thermal radiation in the X-ray 
(below 10.0 keV) and radio domains appears rather unlikely. Our investigation of our XMM-Newton data of NCC 6604 reveals 
a rather sparse distribution of X-ray emitters. Including the two bright non-thermal radio emitters HD 168112 and HD 167971, 
we present a list of 31 X-ray sources along with the results of the cross-correlation with optical and infrared catalogues. A 
more complete spectral analysis is presented for the brightest X-ray sources. Some of the members of NCC 6604 present some 
characteristics suggesting they may be pre-main sequence star candidates. 
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1. Introduction 

The nature of the X-ray emission from hot stars has been the 
subject of many discussions over the last three decades. The 
bulk of the X-ray emission of single O-stars arises in a hot op¬ 
tically thin thermal plasma inside the stellar wind. These plas¬ 
mas are believed to be heated up to several 10® K by shocks re¬ 
sulting from instabilites related to the line driving mechanism 
responsible for the mass loss of massive stars (see e.g. Owocki 
& Rvbicki [19851 Feldmeier et al. ll997> . The case of massive 
binaries is more complex, because they harbour a wind interac¬ 
tion zone that can produce a substantial additional X-ray emis¬ 
sion (Stevens et al. ll992t . Such a wind collision is expected to 
heat the shocked plasma to temperatures of several 10^ K. 


* Based on observations with XMM-Newton, an ESA Science 
Mission with instruments and contributions directly funded by ESA 
Member states and the USA (NASA). Partly based on observations 
collected at the European Southern Observatory (La Silla, Chile). 

** Research Associate FNRS (Belgium) 


Some massive stars are known to display non-thermal emis¬ 
sion in the radio domain (see e.g. Bieging et al. l 19891 Williams 
I1996> . Dougherty & Williams J2000> showed that most of the 
non-thermal radio emitting Wolf-Rayet (WR) stars are bina¬ 
ries, suggesting that the non-thermal phenomenon might be in¬ 
timately related to multiplicity. However, the situation for O- 
stars is much less clear as the binary fraction among O-type 
non-thermal radio emitters is apparently lower than for WR 
stars (see Ranw 120041 . 

The non-thermal radio emission, which is thought to be 
synchrotron emission (White I1985> . requires the presence of 
a population of relativistic electrons inside the radio emitting 
region, as well as the existence of a moderate magnetic field. 
On the one hand, high energy electrons could be accelerated 
through the first order Fermi mechanism in shocks (Pollock 
[TMtI Chen & White IT^ Eichler & Usov [1993}. On the 
other hand, the existence of magnetic fields of a few hundred 
Gauss has recently been confirmed in the case of a few 
massive stars (Donati et al. l2001l I2002t . For a discussion 
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of the physical processes involved in this scenario, see e.g. 
De Becker et al. J2005t . If the ingredients required for non- 
thermal radio emission are present, one may wonder whether 
the relativistic electrons could produce a signature at other, 
essentially higher, energies. Indeed, the intense UV flux from 
massive stars could interact with relativistic electrons through 
Inverse Compton (IC) scattering resulting in a power law 
X-ray emission component (Eichler & Usov 119931 Chen & 
White [T^ . Such a scenario could apply both to single and 
binary systems, since both types of objects are believed to 
harbour hydrodynamic shocks that could accelerate particles. 
There are two major questions to be addressed in order to 
fully understand non-thermal phenomena in early-type stars: 
(1) one has to establish whether non-thermal X-ray emission 
really occurs, and ( 2 ) one has to check whether single stars are 
indeed able to produce non-thermal radiation. 

In the framework of these considerations, non-thermal 
radio emitting massive stars are a priori privileged targets 
to search for an X-ray counterpart to this non-thermal radio 
emission. We therefore initiated a joint X-ray and radio 
campaign aiming at a better understanding of the emission 
processes in non-thermal radio emitters. The first target of this 
campaign observed with XMM-Newton was the 04V((f''')) star 
9 Sgr (Rauw et al. l2002> . Whilst the X-ray spectrum of 9 Sgr 
displays indeed a hard emission tail, the very nature of the 
latter, thermal or non-thermal, could not be fully established. 
Optical spectra of 9 Sgr revealed long term radial velocity 
variations suggesting that 9 Sgr could indeed be a long-period 
binary system. The second object considered was HD 168112 
(De Becker et al. l2004bt . This 05.5III(f+) star, belonging to 
the NGC 6604 cluster (Barbon et al. l2000t . also shows some 
evidence for binarity although no radial velocity variations 
were found. Indeed, De Becker et al. J2004bt report changes 
of the nature of the radio emission (going from composite to 
strongly non-thermal, with a flux variation by a factor 5-7) 
simultaneously with a significant decrease of the X-ray flux 
(about 30%), suggesting that a wind collision occurs in an 
eccentric long period binary system seen under low inclination. 
Blomme et al. ( 12005 at suggest a ~ 1.4yr period for the radio 
fluxes. These results lend further support to the idea that 
binarity might be a necessary condition for non-thermal radio 
emission, also in O-stars. This paper is mainly devoted to the 
study of HD 167971, another non-thermal radio emitter which 
is known to be a multiple system (Leitherer et al. ll987t . and 
belongs also to NGC 6604. 

HD 167971 was first classified as an 08f star (Hiltner 
119561 . This system has been a target of several photometric 
studies (i.e. Johnson [l965l Moffat & Vogt 119751 Leitherer & 
Wolf 119^ . The IR study of Bertout et al. ( I1985> yielded 
M ~ 9.4 X lO“®M 0 yr“^ and Voo ~ 2120kms“^. In a 
more recent photometric and spectrometric study, Leitherer et 
al. ( 119871 proposed this system to consist of a close eclipsing 
binary made up of 2 similar O-stars, with a third more lumi¬ 
nous and more distant O-type companion (05-8V + 05-8V + 
(081)). Their analysis of the light curve revealed a 3.3213 days 
period for the eclipsing binary. They found a mass loss rate of 


~2 X lO“®M 0 yr“^ and a terminal velocity of 3100kms“^ 
for the brightest third star. The study of Davidge & Lorbes 
( I1988t confirmed the orbital period of the close binary sys¬ 
tem. Their photometric solution, which is in better agreement 
with observations, has non-zero third light, suggesting that 
HD 167971 is indeed a multiple system. However, we note that 
at this stage there is no evidence of a gravitational link between 
the third star and the eclipsing binary. 

In his review of radio emitting hot stars, Williams (I1996t 
lists HD 167971 as a non-thermal emitter, with a radio spectral 
index of —0.6 and a 6 cm flux level ranging between 13.8 
and 17.0 mJy. It is the strongest non-thermal radio emitter 
reported in the study of Bieging et al. ( I1989> . This star appears 
consequently to be a good candidate to study the implication 
of binarity on non-thermal emission processes. This paper 
is devoted to the analysis of two XMM-Newton pointings 
obtained in 2002. The radio results from the VLA, as well as a 
discussion of archive radio data, will be presented in a separate 
paper (Blomme et al. l2005bl in preparation). 

NGC 6604 is a rather compact open cluster lying at the 
core of the HII region S54 (Georgelin et al. ll973> . in the 
SerOB2 association (Lorbes l2000t . The study of Lorbes & 
DuPuy ( I1978I I revealed that it is rather young (^4Myr) as 
suggested also by the possible presence of pre-main sequence 
(PMS) objects. However, the study of Barbon et al. (I20()0t did 
not confirm the existence of pre-main sequence stars, although 
these authors also derived an age of the cluster of about 
5 ± 2 Myr. The same authors derive a distance to the cluster of 
about 1.7kpc. In addition to the study of the two non-thermal 
radio emitters mentioned hereabove, the X-ray observation 
of NGC 6604 offers also the possibility to study the stellar 
population of the cluster. Indeed, recent XMM-Newton obser¬ 
vations of NGC 6383 (Rauw et al. lMut and NGC 6231 (Sana 
et al. l2005t reveal a strong concentration of X-ray selected 
PMS stars around the most massive cluster members that are 
found in the cluster core. Such mass segregated clusters with 
a concentration of low-mass PMS objects could indicate that 
the most massive stars have formed through a combination of 
gas accretion and stellar collisions in the very dense core of 
the cluster (Bonnell et al. ll998l Bonnell & Bate l2002l Bonnell 
et al. l2003t . The XMM-Newton observation of NGC 6604 is 
expected to bring new insight on the existence of low-mass 
PMS stars in this cluster. 

This paper is organized as follows. In Sect.|2l we briefly 
discuss the data reduction. Section|3describes the X-ray spec¬ 
tral analysis, as well as some considerations relevant to the 
colliding winds context. In Sect.l^we compare our results for 
HD 167971 with archival X-ray observations. Section|3is de¬ 
voted to a general discussion of HD 167971. The X-ray emis¬ 
sion from the other sources in the held of view is discussed in 
Sect.|^ Linally, our conclusions are presented in Sect.0 

2. Observations 

We obtained two XMM-Newton (Jansen et al. l2001t obser¬ 
vations of the NGC 6604 open cluster during revolution 


























































M. De Becker et al.: An XMM-Newton observation of NGC 6604 


3 


426, in April 2002 (Obs.ID0008820301, JD2452372.477 
- 2452 372.637), and revolution 504, in September 2002 
(Obs.ID0008820601, JD2452526.694 - 2452526.868). The 
aim-point of these observations was set in such a way as to 
monitor simultaneously the two non-thermal radio emitters 
HD 168112 and HD 167971 within a single EPIC field of view. 
The exposure time for both observations was about 13 ks. For 
details on these exposures and on the reduction procedure, 
we refer to the paper devoted to HD 168112 (De Becker et 
al. l2004b> . 

The X-ray events of HD 167971 were selected from inside 
a 60 arcsec radius circular region centered on the star. The 
background was extracted from an annular region around the 
source. The inner and outer radius of this annulus were chosen 
to obtain a surface area roughly equal to the source region area. 
Both regions are shown in Fig.[T]in the case of the three EPIC 
fields, for both observations. In each case, the source region is 
affected by at least one gap, which was excluded by means of 
a rectangular box. These boxes were adjusted following a care¬ 
ful inspection of the adequate exposure maps. Unfortunately 
these gaps sometimes cross the source region close to its cen¬ 
ter, strongly affecting the results discussed hereafter. In the case 
of the EPIC-pn data sets, the background region is crossed by 
a bad column which was also removed with a box. 

3. The EPIC spectrum of HD 167971 

3.1. Spectral analysis 

As already stated by De Becker et al. J2004bt . the end of the 
April 2002 exposure is affected by a strong soft proton flare 
(see Lumb l200^ . However, it appears that the background sub¬ 
traction provides an efficient correction, and the contaminated 
time interval was therefore not rejected. Indeed, the level of the 
background is at least a factor 2 lower than the mean count rate 
inside the source region^. 

To carry out this spectral analysis, we followed the same 
procedure as already applied to other targets like HD 159176 
(De Becker et al. l2004al l and HD 168112 (De Becker et 
al. l2004bt . We chose a list of models including thermal and 
non-thermal emission components. Massive stars are known to 
display thermal spectra with strong and broad emission lines 
(e.g. C Pup, Kahn et al. l20011 t Sco, Cohen et al. l2003t . Most 
current models consider that, in the case of single O-stars, such 
hot plasmas are produced by strong shocks inside the stellar 
wind (Feldmeier et al. ll997t . To model the emission spectrum 
arising from such plasmas, we used an optically thin thermal 
plasma mekal model (Mewe et al. l 19851 Kaastra ll992t avail¬ 
able within XSPEC. In the case of massive binaries, colliding 
winds can also heat the plasma in the wind interaction zone to 
a temperature of up to a few times 10^ K. The emission from 
the wind interaction zone can also be modelled with optically 
thin thermal plasma codes (see e.g. the case of HD 159176, De 
Becker et al. l2004al l. Moreover, to account for the possibility of 

* In the case of HD 168112 (De Becker et al. l2004bV the soft pro¬ 
ton flare did not affect significantly the spectral analysis although the 
background level was of the same order of magnitude as the source 
count rate. 


non-thermal emission arising from the Inverse Compton (IC) 
scattering of UV photons due to the presence of a population 
of relativistic electrons (Pollock 119871 Chen & White 119911 
Eichler & Usov ll993t . power law models were also considered. 
Our grid of models therefore includes pure thermal models, and 
thermal plus non-thermal models. 

The minimization technique was used to estimate the 
quality of the fits described in the following paragraphs of 
this section. The parameter values quoted in Tables |2] and 
0 are obtained through this method. However, to check the 
consistency of our approach, we confronted these results with 
those obtained with other techniques. Indeed, when one deals 
with small count numbers the parameter estimation could be 
achieved through the C-statistic (Cash I1979> . or alternatively 
with the statistic using a Churazov weighting (Churazov et 
al. 119961 . Both methods yield results compatible with those 
quoted in Tables |2]and|3 within the error bars. 


3.1.1. Wind absorption 

In each case, two separate absorption columns were used to ac¬ 
count for interstellar (ISM) and circumstellar (wind material) 
absorption respectively. The first absorption component is due 
to neutral material and the hydrogen column density was fixed 
at 0.63 X 10^^ cm“^. This value was obtained through the gas 
to dust ratio given by Bohlin et al. lITWSl . using a {B — V) 
color index of h- 0.77 (Chlebowski et al. ll989> and an intrinsic 
color index {B — V)o of -0.31 estimated for an 08 star (see 
e.g. Mihalas & Binnev n'981t . For the second absorption col¬ 
umn, left as a free parameter while fitting the spectra, we used 
the ionized wind model already presented by Naze et al. J2004t . 
This model considers the absorption by the 10 most abundant 
elements (H, He, C, N, O, Ne, Mg, Si, S and Fe) fixed at their 
solar abundances (Anders & Grevesse ll989t . We used the stel¬ 
lar fluxes from the Kurucz library of spectra, and a standard 
wind velocity law with /3 = 0.8. 

The main problem to model the wind properties of 
HD 167971 with this model is that it only accounts for the 
geometry of a single star wind, although we know that our 
target is a triple system. The choice of the parameters ad¬ 
equate for our purpose is therefore not straightforward. To 
check the sensitivity of the wind opacity calculated by this 
model regarding the stellar parameters, we computed sev¬ 
eral opacity tables covering a parameter space typical for a 
range of spectral types including that of the components of 
HD 167971. This grid was calculated for mass loss rates rang¬ 
ing from 10 “® to lO“®M 0 yr“^, and terminal velocities be¬ 
tween 2200 and 3100km s“^. These domains cover differ¬ 
ent values proposed by some authors for the parameters of 
HD 167971 (Bertoutet al. [T^ Feitherer et al. [T^ . We used 
also three different Kurucz spectra, respectively referred to as 
low {Teg=30 000; log g=3.5), medium (Tefi=35 000; log 5 = 4 . 0 ), 
and high (reff=40 000; log 5 = 4 . 5 ), and we computed the optical 
depth (t) as a function of the energy for X-ray emitting shells 
located between 1.5 and 100 stellar radii. The opacities result¬ 
ing from the normalization of t by the hydrogen column of the 
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See GIF figure (figl.gif) 


Fig. 1. Source (circle) and background (annulus) regions selected for the spectrum extraction of HD 167971. Boxes were used 
to exclude the CCD gaps and bad columns. We see clearly that some observations are strongly affected by the presence of these 
gaps. The images are displayed in detector coordinates. 


Table 1. Parameters selected to compute the opacity table used 
for the modelling of the absorption of X-rays by the ionized 
wind material. 


Teff (K) 

35 000 

logs 

4.0 

M (Moyr"^) 

5 X 10"® 

Voo (kms“^) 

3100 

Rshell (R*) 

5 


wind is very similar whatever the case considered within this 
parameter space. Consequently, any set of parameters located 
within the boundaries discussed hereabove appears to be ade¬ 
quate. In the next section, we will discuss spectral fittings us¬ 
ing the opacities we obtained from the high Kurucz spectrum, 
and the parameters quoted in Table[2 The parameters were se¬ 
lected on the basis of the values used or derived by Leitherer et 
al. l l 19871 in the computation of synthetic Ha profiles for this 
star. Fig.|2]shows the normalized r as a function of the energy 
for the ionized wind (solid line), along with the same quantity 
for neutral material (dashed line). We see that the largest devi¬ 
ation between the neutral and ionized cases is observed at low 
energies. The opacity table obtained with the parameters listed 
in Tabled was converted into a FITS table in a format suit¬ 
able to be used as a multiplicative absorption model within the 
XSPEC software. It will be referred to as the wind absorption 
model in the remainder of this paper. 


10“®‘ 

“g 10-"" 

O 
3= 

IZ, 

V IQ-"" 

IQ-"^ 

0.5 ] 5 10 

Energy (keV) 

Fig. 2. Opacity resulting from the normalization of the opti¬ 
cal depth by the hydrogen column density as a function of the 
energy for ionized (solid line) and neutral (dashed line) wind 
material. 


3.1.2. Spectral fittings 

As previously mentioned, most of our data sets are strongly af¬ 
fected by the presence of gaps and/or bad columns. Only the 
EPIC-MOSl and EPIC-MOS2 observations of April, and to 
some extent the EPIC-pn observation of September appear to 
be only slightly affected. In all other cases, the CCD gaps lie 
very close to the center of the PSP of the source. Since the data 
seriously affected by gaps yield unreliable spectral fits, we con¬ 
sider only the spectral fits for those data that are least affected. 

Two-component models were used as a first step with some 
success, but it appeared that those models were unable to fit 
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Table 2. Parameters for EPIC spectra of HD 167971 in the case of a wabs*wind* (mekal+mekal+mekal) model. Results 
are given for MOSl, MOS2, combined MOS in the case of the April observation, and for EPIC-pn only in the case of the 
September one. The hrst absorption component is frozen at 0.63 x 10^^ cm“^. The second absorption column, quoted as (in 
cm“^), stands for the absorption by the ionized wind material. The last column gives the observed flux between 0.4 and 10.0 
keV. The normalization parameter (Norm) of the mekal component is dehned as (10“^^/(4 tt D^)) f Ug nn dV, where D, rig 
and riH are respectively the distance to the source (in cm), and the electron and hydrogen number densities (in cm“^). The error 
bars represent the l-u conhdence interval. 



Log Aw 

kTi 

(keV) 

Normi 

feTa 

(keV) 

Norma 

kTi 

(keV) 

Norms 

xl 

d.o.f. 

Obs.Flux 
(ergcm-2 s-*) 

April 

MOSl 

21.77 

±0.08 

0.27 

±0.02 

1.38 X 10~^ 

± 0.67 X 10"^ 

0.79 

±0.08 

2.05 X 10-® 

± 0.72 X 10-® 

3.30 
± 1.32 

6.38 X 10-“^ 

± 3.86 X 10-* 

0.867 

120 

1.79 X 10-*2 

MOS2 

21.78 

±0.07 

0.24 

±0.02 

1.68 X 10“^ 

± 0.76 X 10-2 

0.77 

±0.06 

2.62 X 10-® 

± 0.76 X 10-® 

3.96 

±2.68 

4.26 X 10-* 

± 3.81 X 10-* 

1.115 

116 

1.69 X 10-*2 

MOSl 

-I-MOS2 

21.78 

±0.05 

0.25 

±0.01 

1.52 X 10-2 
± 0.50 X 10-2 

0.78 

±0.05 

2.35 X 10-® 

± 0.52 X 10-® 

3.57 
± 1.25 

5.21 X 10-* 

± 2.73 X 10-* 

0.973 

243 

1.70 X 10-*2 

September 

pn 

21.69 

±0.06 

0.27 

±0.01 

9.55 X 10-® 

± 2.86 X 10-® 

0.77 

± 0.05 

1.87 X 10-® 

± 0.43 X 10-® 

2.40 

±0.43 

5.47 X 10-* 

± 2.32 X 10-* 

1.089 

244 

1.49 X 10-*2 


Table 3. Same as Table|2] but in the case of a wabs*wind* {mekal+mekal+power) model. Eor the power law component, 
the normalization parameter (Norma) corresponds to the photon flux at 1 keV. 



Log Aw 

kTi 

(keV) 

Normi 

fcTa 

(keV) 

Norma 

r 

Norms 

xl 

d.o.f. 

Obs.Flux 
(ergcm-^ s-*) 

April 

MOSl 

21.48 

±0.28 

0.31 

±0.05 

4.24 X 10-® 

± 5.77 X 10-® 

0.75 

±0.08 

1.36 X 10-3 
± 0.78 X 10-3 

3.04 
± 0.282 

8.47 X 10-* 

± 5.09 X 10-* 

0.851 

120 

1.70 X 10-*^ 

MOS2 

21.66 

±0.14 

0.25 

±0.02 

9.54 X 10-^ 

± 6.66 X 10-3 

0.75 

±0.06 

2.04 X 10-3 
± 0.94 X 10-3 

2.95 
± 0.426 

6.07 X 10-* 

± 6.25 X 10-* 

1.095 

116 

1.63 X 10-*^ 

MOSl 

-tMOS2 

21.63 

± 0.11 

0.27 

±0.02 

8.02 X 10-3 
± 4.65 X 10-3 

0.76 

±0.05 

1.78 X 10-3 
± 0.64 X 10-3 

2.94 

± 0.268 

6.78 X 10-* 
±4.21 X 10-* 

0.954 

243 

1.70 X 10-*^ 

September 

pn 

21.68 

±0.07 

0.27 

±0.01 

9.08 X 10-3 
±3.18 X 10-3 

0.78 
± 0.05 

1.84 X 10-3 
± 0.47 X 10-3 

2.51 

±0.33 

2.96 X 10-* 

± 2.78 X 10-* 

1.072 

244 

1.52 X 10-*^ 


data above about 4-5 keV. As a consequence, we used three 
component models, each one having two mekal thermal com¬ 
ponents, plus another thermal or power law component. Solar 
abundances are assumed throughout our htting procedure. 
The results obtained with these two models are summarized 
in Tables |2] and 0 The last column of these tables yields 
the observed (i.e. absorbed) flux integrated between 0.4 and 
10.0 keV for each set of model parameters. The dispersion on 
the flux was evaluated from the range of fluxes for the models 
with parameters covering the confidence intervals quoted in 
Tables |2] and 0 According to this approach, the error on the 
flux in that energy band is expected to be about 5 %. 

Eirst, the EPIC-MOS1 and EPIC-MOS2 results of the April 
observation appear to be very consistent, and our models were 
also fitted to the two data sets simultaneously. A comparison 


between the results of the two models reveals no significant dif¬ 
ferences in the quality of the fit. Both models are nearly equiva¬ 
lent within the error bars for the two first thermal mekal com¬ 
ponents. The values obtained and used as goodness-of-fit 
criterion are very similar for all EPIC instruments, either con¬ 
sidered independently or combined. Moreover, the fluxes esti¬ 
mated on the basis of these two models are similar within the 
error bars, reinforcing the idea that both models reproduce the 
spectral shape rather consistently. Unfortunately, the poor qual¬ 
ity of the data above 5 keV does not allow us to unambiguously 
check for the presence of the Ee K line that could help constrain 
the nature of the hard spectral component. These results show 
that the ambiguity of the nature of the hard part of the spec¬ 
trum (thermal versus non-thermal) can not be alleviated with 
our data. This case is reminiscent of the results of the XMM- 
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EPIC-M0S1 ond EPIC-M0S2: 3-T MODEL 



0.5 1 2 5 10 

channel energy (keV) 


Fig. 3. Combined EPIC-MOSl and EPIC-MOS2 spectra of 
HD 167971 for the April 2002 observation, fitted with 
a wabs*wind* (mekal+mekal+mekal) model between 
0.4 and 10.0 keV. 


EPIC —pn; 3—T model 



0.5 1 2 5 10 

channel energy (keV) 


Fig. 4. EPIC-pn spectrum of HD 167971 of 
the September 2002 observation, fitted with a 
wabs*wind* (mekal+mekal+mekal) model between 
0.4 and 10.0 keV. The three components are individualy 
displayed. 


Newton observations obtained in the case of other non-thermal 
radio emitters (see e.g. 9 Sgr, Rauw et al. l2002t . 

The EPIC-pn observation of September gives results very 
similar to those of the EPIC-MOS ones of April. Even if some 
systematic deviations could be expected as the EPIC-MOS and 
EPIC-pn instruments are somewhat different, we see that the 
overall spectral shape remained more or less steady between 
the two observations. If we consider the three-temperature 
model, we see that the emission measure of the soft ther¬ 
mal component decreases from April to September. This de¬ 
crease is responsible for the lower flux quoted in Table|2]for 
the September pointing. However, if we consider the model de¬ 
scribed in Tabled it is the normalization factor of the power 
law which decreases between the two pointings whilst the pa¬ 
rameters of the two thermal components remain steady. 


Table 4. Observed count rates for the three EPIC instruments 
evaluated between 0.4 and lO.OkeV. The error bars on the count 
rates represent the l-tr confidence interval. 


Observation 

MOSl CR 
(cts s“^) 

MOS2 CR 
(cts s“^) 

pnCR 

(ctss“^) 

April 2002 

0.146 

± 0.004 

0.143 

± 0.004 

- 

September 2002 



0.465 
± 0.008 


3.2. X-ray luminosity of HD 167971 
3.2.1. X-ray variability 

When investigating the X-ray variability of a system like 
HD 167971, one has to consider several scenarios likely to 
produce variations in the X-ray luminosity. Variability could 
indeed be expected for several reasons. Eirst, the interaction of 
the stellar winds of the two components of the 05-8V + 05-8V 
close binary system could vary with the phase. Since the close 
binary has a circular orbit (van Genderen et al. ll988t . the 
only way to produce a modulation of the X-ray flux from the 
wind-interaction zone would come from a change in the col¬ 
umn density along the line of sight as a function of the orbital 
phase. The resulting change in absorption would then mostly 
affect the soft part of the EPIC spectrum. Nevertheless, no 
significant change appears in the spectrum or in the absorption 
component of the model used to fit the data between the first 
and the second exposure. A second possible origin could be the 
putative wind-wind interaction due to the presence of the third 
star, if the system in indeed triple. The stellar wind of the third 
component is indeed likely to produce a second interaction 
zone by colliding with the winds of the stars of the close binary 
system. Einally, variability could possibly result from the 
eclipse of the 05-8V + 05-8V binary system. According to the 
ephemeris derived by Leitherer et al. (11987), the April obser¬ 
vation was performed between phases 0.64 and 0.69, and the 
September one between phases 0.04 and 0.12. As phase zero 
corresponds to the primary eclipse, we see that the September 
XMM-Newton observation was performed when the primary 
was partially hidden by the secondary, allowing us to possibly 
expect a lower X-ray flux than during the April observation 
when no eclipse occurred. However, the extrapolation of this 
ephemeris to the epoch of our observations must be considered 
with caution, especially for such a short period. To clarify this 
question, we investigated the X-ray variability of HD 167971 
on the basis of the results of our two XMM-Newton obser¬ 
vations, along with a confrontation to archive data (see Sect.|4}. 

Count rates of HD 167971 between 0.4 and 10.0 keV are 
presented in Table|3 Since CCD gaps (see Eig.Q affect the 
area of the effective source region considered for the event se¬ 
lection, the variability between our two observations can not be 
discussed on the basis of count rates. 

To obtain a more consistent, although model dependent 
way to quantify the X-ray emission, absorbed fluxes and fluxes 
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corrected for interstellar absorption were evaluated on the 
basis of the pure thermal model. The fluxes quoted in Table|2] 
give a rough idea of the behaviour of HD 167971 between 
April and September 2002. If we make the assumption that 
we can compare the fluxes obtained separately with the 
EPIC-MOS and EPIC-pn instruments, we observe a slight 
decrease of the X-ray flux (~12%). However, we remind 
that the error on the flux estimate is about 5 %, and that 
such a heterogeneous comparison may introduce another 
systematic error of a few percent. Therefore, the variability 
of HD 167971 between the two XMM-Newton observations 
is of rather low level. If this slight variation is real, it can be 
due to the eclipse occurring during the second observation as 
pointed out above. We note however that the decrease of the 
emission measure is observed mostly for the soft component 
in the case of the three-temperature model, whilst the case of 
V444Cyg discussed by Pittard J2002t shows that the eclipse 
would indeed be expected to reduce mainly the flux in the 
hard part of the X-ray spectrum. A detailed modelling of the 
effect of the absorption by the winds, of the occultation of 
the colliding zone, and of the dynamic of the stellar winds 
including radiative inhibitaion effects is needed to investigate 
the effect of the eclipse on the observed X-ray spectrum. 

The variability on shorter time scales was also investigated. 
Light curves were generated in different energy bands between 
0.4 and lO.OkeV, with time bins ranging from 100 s to 1000s. 
These light curves were background corrected, and each time 
bin accounts for Good Time Intervals (GTIs)^. Variability tests 
(chi-square, /7ov-test^) were applied to all of them and no sig¬ 
nificant short term variability was found. 

3.2.2. Overall luminosity 

To obtain an estimate of the overall X-ray luminosity corrected 
for the ISM absorption, let us consider the model parameters 
obtained in the case of EPIC-MOS for the first observation, 
and EPIC-pn for the second one. This yields an Lx of about 
4.3 X 10^^ and 3.8 x 10^^ ergs“^ respectively for the April and 
September observations, for a distance of 2 kpc. According to 
the bolometric magnitudes given by van Genderen et al. J1988t 
for the three components of this presumably triple system, 
we can infer bolometric luminosities of 1.5 x 10^® ergs“^ for 
each star of the close binary system, and 3.0 x 10^® erg s“^ for 
the third component. This yields a total bolometric luminos¬ 
ity of about 6.0 X 10^®ergs“^. Eollowing the empirical rela¬ 
tion given by Berghofer et al. \\991\ applied to the three mem¬ 
bers of the multiple system, we infer expected individual X- 
ray luminosities and add them to obtain a total X-ray luminos¬ 
ity of about 1.0 X 10^^ergs“^ corresponding to an expected 
LxILhoi ratio of about 1.7 x 10“^. Eollowing the X-ray lumi¬ 
nosities derived hereabove, the observed Lx/Lboi ratio goes 

^ Even if no flare contaminated time interval was rejected, standard 
GTIs are anyway always applied to XMM-Newton data, and should be 
taken into account in every timing analysis. 

^ See Sana et al. <2004> for a detailed discussion of this variability 
test. 


0 ) 



Date 


Fig. 5. Normalized equivalent X-ray count rates of HD 167971 
arising from different observatories, as a function of time. 
1; EINSTEIN-\VC, March 1981. 2: ROSAT-VSVC, September 
1993. 3: /?6)5Ar-HRI, September and October 1995. 4: XMM- 
Newton-¥P\C, April 2002. 5: XMM-Newton-¥P\C, September 
2002. We note that the ROSAT and EINSTEIN points are aver¬ 
aged values from several observations spread over a few hours 
to a few days. 


from 7.2 x 10“^ to 6.3 x 10“^ respectively for the April and 
the September XMM-Newton observations. These values yield 
X-ray luminosity excesses of about a factor 4.2 and 3.7 for both 
observations respectively. These values are significant consid¬ 
ering the scatter of the relation of Berghofer et al. M9911 . sug¬ 
gesting that HD 167971 is overluminous in X-rays. 

4. Comparison with previous observations 

HD 167971 was observed with the EINSTEIN satellite on 31 
March 1981 (sequence number 4240, 5.22ks). Chlebowski et 
al. if!^ reported an X-ray luminosity of 5.8 x 10^^ ergs ^ 
for a 2 kpc distance. This value is slightly larger than owe XMM- 
Newton results. However, we compared the EINSTEIN results 
to ours following a more model independent way. We passed 
the model obtained in the case of the three thermal mekal 
components through the LLV^LELV response matrix to obtain a 
simulated IPC spectrum. The model parameters are taken from 
the simultaneous EPIC-MOS fit for the April observation, and 
from EPIC-pn for the September one (see TableEJi. Count rates 
were estimated between 0.2 and 3.5 keV (the energy range of 
the IPC instrument) on the basis of these simulated spectra, and 
we obtained 0.049 ± 0.001 ctss“^ and 0.047 ± 0.001 ctss“^ 
respectively for the April and the September 2002 exposures. 
These values are in excellent agreement with the count rate 
given by Chlebowski et al. lfT98^ . i.e. 0.050 ± 0.005ctss-^. 
This result suggests that at the time of the EINSTEIN observa¬ 
tion, HD 167971 was in an emission state similar to owe XMM- 
Newton observations. 

This target was also observed with ROSAT. A PSPC 
observation (rp500298n00, 9.29 ks, performed between 13 
September and 15 September 1993) reports a count rate of 
0.129 ± 0.004cts s“^. By folding the three-temperature model 
obtained for the fit of our EPIC data through the PSPC re¬ 
sponse matrix, we obtain 0.106 ± 0.002ctss“^ and 0.103 ± 
0.002 cts s“^ respectively for our two XMM-Newton observa¬ 
tions between 0.1 and 2.5 keV (the PSPC bandpass). These re- 
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suits suggest that HD 167971 was in a higher emission state at 
the time of the ROSAT observation. The same procedure was 
applied in the case of a HRl observation performed between 
12 September and 10 October 1995 (rh201995n00, 36.3 ks). 
The reported HRl count rate is 0.035 ± 0.001 cts s“^. The re¬ 
sults we obtain by passing the model through the HRl response 
matrix are 0.042 ± 0.001 cts s“^ and 0.041 ± 0.001 cts s“^ re¬ 
spectively for the April 2002 and the September 2002 XMM- 
Newton observations, suggesting that at the moment of the HRl 
observation, HD 167971 was in an emission state somewhat 
lower than that observed with XMM-Newton. 

In summary, these results from archive data and from our 
XMM-Newton observations are plotted in Fig.|3 We note that 
we do not expect any confusion with other X-ray sources to 
affect the archive count rates quoted in this section. The count 
rates are normalized relative to the XMM-Newton count rate of 
April 2002 arbitrarily set to unity. The normalized count rate of 
the September 2002 XMM-Newton observation is determined 
following the flux ratio observed between 0.4 and 10.0 keV 
for the three-temperature model. The error bars on the XMM- 
Newton count rates are set to about 5 % in Fig.|4] We consider 
that this value is more realistic than the standard deviations 
quoted in Tabled We see that the X-ray flux of HD 167971 un¬ 
dergoes a significant variability with a peak to peak amplitude 
of about 40 %. However, we note that this apparent variabil¬ 
ity relies mostly on the position of the ROSAT-Y’SA’C point in 
Fig.|5l as other count rates suggest a more or less constant level. 
Moreover, we want to draw attention to the fact that the overall 
shape of the light curve plotted in Fig.|3is very similar to that 
presented in Fig. 9 of De Becker et al. J2004bt for HD 168112, 
even if the amplitude of the apparent variability was larger in 
that case. Consequently this apparent variability might be at¬ 
tributed to some systematic effects affecting the various obser¬ 
vations (common to HD 168112 and HD 167971) mentioned in 
this section. A detailed description of this behaviour requires a 
better data coverage than that available at the time of this study. 

5. Discussion 

5.1. Search for a colliding-wind signature 

As HD 167971 is a multiple system whose components are 
O-stars, a wind-wind interaction is expected. Do our XMM- 
Newton data reveal the signature of colliding winds? Several 
questions should be considered to address this point; 

- Plasma temperature. The results of the fittings with the 
three-temperature model point to a plasma temperature 
for the hardest component of about 2 - 4keV, i.e. about 
2.3 X 10^ to 4.6 X 10^ K. Such high temperatures are usu¬ 
ally not expected for the emission from shocks due to in¬ 
trinsic wind instabilities, but are expected for shocks be¬ 
tween the winds within a binary system (see e.g. Stevens 
et al.[T993. These post-shock temperatures correspond to 
pre-shock velocities of about 1300 - 1800 kms“^, typ¬ 
ical for winds colliding at speeds near the terminal ve¬ 
locity, as should be expected for long period wide bina¬ 
ries (Pittard & Stevens fl997> . For short period binary sys¬ 
tems like the 05-8V H- 05-8V pair in HD 167971 (P ~ 


3.32 d), the winds have not reached their terminal ve¬ 
locities before they collide, and the temperatures are not 
expected to be so high. Such systems yield typical tem¬ 
peratures of at most 1 X 10^ K. For instance, HD 159176 
(07V H- 07V) with a period of 3.367 d displays tem¬ 
peratures of about 2 X 10® to 6 x 10®K (De Becker et 
al. l2004at . whilst the 07.5(f)III H- 07.5(f)III close binary 
HD 152248 (5.816 d period) has plasma temperatures lower 
than 1 X 10^ K (Sana et al. l2004t . So, if the high tempera¬ 
tures we derived from our spectral fittings have a physi¬ 
cal meaning, they are most probably related to the interac¬ 
tion between the wind of the 081 component and the com¬ 
bined winds of the close 05-8V H- 05-8V binary system. 
We note however that high plasma temperatures can arise 
in the winds of individual stars if the magnetic field of the 
star is able to confine the stellar wind near the magnetic 
equatorial plane (Babel & Montmerle ll997> . 

- X-ray luminosity. The luminosities we derived from our 
spectral fittings point to an excess attributable to a wind- 
wind interaction. This is compatible with a scenario where 
the colliding winds bring a significant contribution to the 
overall X-ray luminosity, in addition to the emission from 
the intrinsic shocks of the individual winds of the compo¬ 
nents of HD 167971. 

- X-ray variability. As discussed in Section l3.2.1l the ques¬ 
tion of the variability of the X-ray flux is rather complex. 
We do not detect any strong variability in the X-ray flux 
between our two XMM-Newton observations. Only a slight 
decrease possibly due to the eclipse of the 05-8V ■¥ 05- 
8V close binary system is observed. However, our discus¬ 
sion of archive data reveals that HD 167971 possibly shows 
some variability, even if the data coverage is insufficient 
to constrain the time scale of this variability. If not due to 
some undetermined systematic effect (see last paragraph of 
Sect.|3l, the amplitude of this apparent variability seems to 
be rather high (~ 40 %). Even if there is no clear evidence 
that the three stars are physically related, this variability 
might be compatible with what could be expected from a 
long period and eccentric binary system. The count rates 
from ROSAT and EINSTEIN result indeed from a combi¬ 
nation of several obervations spanning different phases of 
the orbit of the 05-8V H- 05-8V binary system, and conse¬ 
quently are phase averaged values. Moreover, most of the 
count rates quoted in Fig.|3 point to a low level, although 
it is not reasonable to think that all the pointings except 
that of RD^AT-PSPC fall at the moment of an eclipse. For 
these reasons, we do not expect the eclipse of the close bi¬ 
nary system to be responsible for the variability illustrated 
in Fig.|3 However, considering the lack of reliability of this 
putative long term variability, we decided to attribute only 
little weight to this argument. 

Regarding these results, we realize that we find a probable 
signature of a wind-wind collision. Some clues, like the high 
plasma temperature in the context of the thermal model, point 
to a possible detection of the interaction between the close sys¬ 
tem and the third more luminous companion. With the large 
separations characterizing this interaction, the shocks of the 
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collision zone are expected to be strongly adiabatic, resulting in 
a phase-locked X-ray variability scaling with 1/D, with D be¬ 
ing the distance from one of the two components to the inter¬ 
action zone (see Stevens et al. l 19921 Pittard & Stevens [T997> . 
According to this scenario, the X-ray emission should peak 
close to the periastron passage. Such a high temperature for 
the hard component of the spectrum, if it is thermal, was also 
reported in the case of HD 168112 which is possibly a wide 
and eccentric binary system (De Becker et al. l2004bl Blomme 
et al. l2005at . In the case of HD 167971, the X-ray luminosity 
excess attributed to the wind collision is much more significant. 

Quantifying the respective contributions from the two ex¬ 
pected collision zones is not an easy task. To address such 
an issue, an accurate knowledge of the stellar properties of 
the three stars along with a better idea of their relative posi¬ 
tions is needed. We can however tentatively adopt the following 
semi-quantitative approach. From Sect. 13.2.^ we can assume 
that the cumulated intrinsic contributions from the three stars 
is about 10^^ erg s“^. As the overluminosity factor is about 4, 
the X-rays arising from the collision zone(s) amount to about 
3 X 10^^ ergs“^. On the other hand, for the 07V + 07V close 
binary system HD 159176 (De Becker et al. l2004a^ rather simi¬ 
lar to the eclipsing binary harboured by HD 167971'*, the same 
considerations lead to an X-ray luminosity for the collision 
zone of about 1.3 x 10^^ergs“*. If the 05-8V + 05-8V bi¬ 
nary within HD 167971 produces the same amount of X-rays, 
we estimate that the putative collision zone due to the third 
companion may produce more than twice as much X-rays as 
that of the close eclipsing binary. We mention however that the 
highest temperature component of the thermal model used to fit 
the EPIC spectra accounts for only about 10 % of this quantity, 
suggesting that the softest components account significantly for 
the putative wind-wind collision with the third star. 

5.2. Non-thermal emission from HD 167971 

In the radio domain, the non-thermal nature of the emis¬ 
sion from HD 167971 is well established (Bieging et al. l 19891 
Blomme et al. l2005bt . The high non-thermal radio flux sug¬ 
gests that it comes from the interaction between the close 
eclipsing binary and the third star. However, in the X-ray do¬ 
main, the nature of the hard component of the spectrum is still 
unclear. As discussed in Sect. l3.ll the pure thermal model and 
that including the power law yield results of similar quality. 
This ambiguity is a common feature in the study of X-ray spec¬ 
tra of massive stars (see e.g. De Becker l200H Rauw et al . 120021 
De Becker et al. l2004bll . If we make the assumption that a non- 
thermal emission component is responsible for the hard part 
of the spectrum of HD 167971, we obtain photon index values 
close to 3, which is rather steep as compared to the value, i.e. 
1.5, considered by Chen & White J1991t for X-rays produced 
through Inverse Compton scattering from relativistic electrons 
accelerated in strong shocks. For a discussion of such high val¬ 
ues of the photon index, we refer to De Becker et al. (I2004bt . 

* The bolometric luminosity provided by van Genderen et al. <1988> 
for the latter stars is intermediate between those given by Howarth & 
Prinia fT98^ for the 05.5V and 06V spectral types. 


5.3. Comparison with other eariy-type stars 

Up to now, our campaign devoted to non-thermal radio emit¬ 
ters includes the observations of three targets: 9 Sgr (Rauw et 
al.|2003, HD 168112 (De Becker et al. l2004b> . and HD 167971 
(this study). The next targets of this campaign are the non- 
thermal radio emitters of the Cyg OB2 association (#8 A, #9 and 
#5, see e.g. Waldron et al. ll998t . Forthcoming studies based on 
XMM-Newton and INTEGRAL observations of these stars will 
bring new elements to the overall discussion of the high energy 
emission from massive stars. For instance, CygOB2#8A is a 
well known non-thermal radio emitter that has recently been 
identified as a binary system by De Becker et al. J2004ct . lend¬ 
ing further support to the scenario where binarity is a neces¬ 
sary condition to observe non-thermal emission from massive 
stars. This idea is also supported by the results of Benaglia & 
Koribalski J2005t who discuss the case of four southern non- 
thermal radio emitters, among which three are confirmed bi¬ 
nary systems. The need for a binary scenario to explain the non- 
thermal radio emission has independently been demonstrated 
by Van Loo following a theoretical approach. Several 

works has recently been devoted to the non-thermal radio emis¬ 
sion from massive binaries (Dougherty et al. l2003l Pittard et 
al . 12(1051 . These recent models take into account several phys¬ 
ical effects likely to affect the observed synchrotron spectrum 
(Razin effect, free-free absorption, IC cooling). Future devel¬ 
opments of these models are expected to lead to an estimate of 
the non-thermal high-energy emission from massive binaries. 

At this stage, it could be interesting to confront the results 
obtained in the case of the first three targets studied in this cam¬ 
paign. The summary of the results arising from these studies is 
provided in Tabled For these three objects, our analyses did 
not reveal any unambiguous power-law emission component in 
the X-ray spectra. A striking characteristic of their X-ray spec¬ 
tra is that the fit with thermal models reveals a rather hard emis¬ 
sion, with plasma temperatures of a few 10*^ K. HD 167971 is 
a known multiple system and such high plasma temperature 
could be explained by a wind-wind interaction. HD 168112 is 
not known to be a binary system, but was proposed to be a 
binary candidate by De Becker et al. J2004b> . The results of 
the study of Blomme et al. J2005at lend further support to this 
idea, suggesting a period of about 1.4 yr based on archival ra¬ 
dio and X-ray data, even if no unambiguous evidence of bi¬ 
narity has been found in the optical domain. In the case of 
9 Sgr, low amplitude radial velocity variations in the optical 
spectrum suggest it is a binary with a yet undetermined period 
(Rauw et al. l2002t . If the high plasma temperatures observed 
for these three stars are indeed due to wind collisions within 
binary systems, the observed post-shock temperatures are re¬ 
lated to rather high pre-shock wind velocities (typically at least 
~ 1000km s“*). Such high velocities can only be achieved in 
binary systems where the winds collide after they have reached 
their terminal velocity, i.e. in systems with periods exceeding 
several days. In this scenario, the non-thermal radio emission 
is produced by a population of relativistic electrons accelerated 
in the shock due to the same wind-wind interaction. The pu¬ 
tative non-thermal X-ray counterpart to this non-thermal radio 
emission has little chance to be detected, as it would probably 
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Table 5. Comparison of the cases of 9 Sgr (Rauw et al.|2004j, HD 168112 (De Becker et al. l2004bt . HD 167971 (this study), and 
HD 159176 (De Becker et al. l2004at . The question marks for the multiplicity and the period stem from the fact that these are not 
well established. 



9 Sgr 

HD 168112 

HD 167971 

HD 159176 

Multiplicity 

binary ? 

binary ? 

triple 

binary 

Period 

long period ? 

long period (1.4yr) ? 

3.3213 d 
-1- long period 

3.367 d 

Thermal X-rays 

soft - 1 - hard 

soft - 1 - hard 

soft - 1 - hard 

soft 

Hard X-rays 

ambiguous 

most probably 
thermal 

ambiguous 

possible power 
law tail 

Synchrotron radio emission 

yes 

yes 

yes 

no 


be overwhelmed by the hard thermal component produced by 
the wind collision. If on the contrary the period were shorter, 
the pre-shock velocities would be lower and the X-ray thermal 
emission would be softer, allowing consequently the putative 
non-thermal X-ray emission to be detected in the hard part of 
the spectrum. However, the shorter orbital period required to 
possibly observe the non-thermal X-ray emission could cause 
the wind collision zone to be deeply embedded in the stellar 
wind material, thereby leading to a severe absorption of the 
synchrotron radio emission. 

The fact that the conditions for detecting a non-thermal 
emission in the radio domain may be different from those for 
such a detection in the X-rays raises an interesting question. 
We may indeed wonder whether non-thermal radio emitters are 
ideal candidates to search for a non-thermal X-ray counterpart. 
To address this issue, let us consider some results obtained for 
other massive stars recently observed in the X-ray domain. 

In the case of the short period massive binary HD 159176 
(OV7 + OV7), a high energy tail is observed in the X-ray 
spectrum. It can not be fitted with a thermal model^, but can 
be approximated with a power law (De Becker et al. l2004al . 
For that system, the typical temperature of the X-ray emitting 
plasma is sufficiently low to prevent the putative power law 
tail to be overwhelmed by the thermal X-rays from the collid¬ 
ing wind zone. This is due to the fact that the winds collide 
as they have not reached their terminal velocities. The non¬ 
detection of HD 159176 (Bieging et al . 119891 in the radio do¬ 
main is compatible with the rather low flux expected for purely 
thermal radio emitters, and can be explained by the fact that 
any putative synchrotron emitting zone would be deeply em¬ 
bedded within the stellar winds. We mention also the case of 
WR 110 (HD 165688) whose binarity has not yet been estab¬ 
lished. Skinner et al. J2002t reported a soft thermal X-ray spec¬ 
trum, along with a high energy power law tail, whilst its radio 
spectrum is purely thermal. 

Considering these recent results, we propose that the non- 
thermal radio emission could ideally be detected in the case of 
binaries with periods larger than several weeks. For instance, 
the shortest period well-established non-thermal radio emitter 
is CygOB2#8A, with an orbital period of about 21.9 d (De 
Becker et al. l2004ct . On the other hand, short period binaries 

^ Optically thin plasma model (Mewe et al . 119851 Kaas tra il9921 or 
colliding wind model (Antokhin et al. l2004> . 


are not expected to display a non-thermal radio emission be¬ 
cause of the strong absorption by the wind material. However, 
their rather soft thermal X-ray emission could possibly unveil 
a power law component produced by inverse Compton scatter¬ 
ing of UV photons. Moreover, the less diluted UV radiation 
field in the shock region of close binary systems is expected 
to lead to the production of a higher non-thermal X-ray flux, 
favoring its detection. In the context of this scenario, the si¬ 
multaneous observation of non-thermal radiation in the X-ray 
(below lO.OkeV) and radio domains appears rather unlikely. 

6. Other X-ray sources in the field of view 

6.1. Source list 

Beside HD 167971 discussed in the previous sections, and 
HD168112 (De Becker et al. l2004bl l. other fainter X-ray 
sources are observed in the held of the EPIC cameras. To inves¬ 
tigate the X-ray emission from these sources, we processed our 
data with the version 6.0.0 of the Science Analysis Software 
(SAS) to benefit from the latest version of the edetect_chain 
metatask. We Altered the event lists of the April observation to 
reject the end of the exposure contaminated by a soft proton 
flare (see Sect. 3.1 of De Becker et al. l2004bll . thereby reducing 
the effective exposure time to about 9 and 6 ks for EPIC-MOS 
and EPIC-pn respectively. Adding together the two observa¬ 
tions and considering that EPIC-pn is about twice as sensitive 
as EPIC-MOS, we obtain an EPIC-MOS equivalent combined 
exposure time of about 80 ks on NGC 6604. 

We first applied the source detection simultaneously to the 
data from the three EPIC instruments of our two observations 
to improve the detection efficiency. We used three energy bands 
respectively refered to as S (0.5-1.0keV), M(1.0-2.5keV), and 
H (2.5-lO.OkeV). The images in these three energy bands were 
extracted on the basis of the merged event lists of the two ob¬ 
servations, for the three EPIC instruments respectively. We ob¬ 
tained a first source list using a detection likelihood threshold 
of 40. We then inspected each source of the list by eye to reject 
false detections due to background fluctuations and instrumen¬ 
tal artifacts following the same approach as Rauw et al. (I2003> 
and Naze et al. il2?)n4l . Adopting the formalism described by 
Sana et al. (|2005}, we determined the likelihood thresholds ad¬ 
equate for individual datasets. The likelihood thresholds we 
Anally adopted are 40 for the combined EPIC data sets, 40 
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See GIF figure (fig6.gif) 


Fig. 6. Image of the NGC 6604 open cluster from the combined EPIC data sets of the two observations. The X-ray sources listed 
in Table|6lare labelled. The pixel size is 2.5 ". The right ascension is increasing from the right to the left, and the declination is 
increasing from the bottom to the top. The width of the field is about 30'. 


for combined EPIC-MOS, 10 for individual EPIC-MOS, and 
10 for EPIC-pn. After applying the same procedure to data 
sets resulting from individual instruments and/or observations, 
we added a few more sources to the list obtained for the full 
combined data set. We finally obtain a catalogue of 31 X-ray 
sources presented in Tablel^by order of increasing right ascen¬ 
sion. The position of these sources in the EPIC field is shown 
in Pig.|6l All sources are identified following the naming con¬ 
ventions recommended by the XMM SOC and the lAU. 

6.2. Source identification 

The position of the 31 X-ray sources were cross-correlated 
with three catalogues: the Guide Star Catalogue (GSC, version 
2.2)®, the Two Micron All Sky Survey (2MASS, Skrutskie et 
al. [T997t . and the US Naval Observatory (USNO-Bl.O, Monet 
et al. l2003t . We derived the optimal correlation radius follow¬ 
ing the procedure described by leffries et al. (119971 . In this 


® The Guide Star Catalogue-II is a joint project of the Space 
Telescope Institute and the Osservatorio Astronomico di Torino. 


way, we find that a cross-correlation radius of 4.5 arcsec in¬ 
cludes the majority of the true correlations while rejecting most 
of the spurious correlations. With a radius of 4.5 arcsec we ex¬ 
pect to achieve at least 23 true and only about 2 spurious cor¬ 
relations. We also assume that the position of the source on the 
detector has little impact on the correlation radius. We should 
however keep in mind that this assumption is mostly valid close 
to the center of the field for a given instrument, but could pos¬ 
sibly bias our source identification procedure for large off-axis 
angles. 

Among the 31 X-ray sources quoted in Table|^ five sources 
(#12, #13, #14, #16 and #25) have neither a GSC, 2MASS 
nor USNO counterpart. 20 (respectively 23) have at least one 
optical counterpart in the GSC (resp. USNO) catalogue, and 
some sources have possibly two optical (#20, #27 and #29) 
counterparts. No other positive correlations were found with 
the sources included in the photometric survey of Eorbes & 
DuPuv( ll978t within a 10 arcsec radius. A total of 18 EPIC 
sources in the field of view of NGC 6604 have a single 2MASS 
counterpart while another 6 sources have two infrared coun¬ 
terparts. In Pig.|8] we show the JHK colour-colour diagram 
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Table 6. Characteristics of the X-ray sources in the field of view. The 31 sources are sorted by order of increasing right ascension. Source #1 was identihed as an extended 
source, and the 30 next ones as point sources. The count rates corrected for the exposure map are quoted for all three EPIC instruments. The missing values are due to a strong 
deviation from those of the other instruments because of CCD gaps. The hardness ratios are given for EPlC-pn only. The quoted counterparts are located within 4.5" of the 
X-ray sources. Nr is the number of counterparts within the correlation radius, and d is the angular separation between the X-ray source and its nearest counterpart. The error 
bars on the count rate represent the ± 1 cr Poissonian standard deviation. 


# 

(1) 

XMMU 

(2) 

MOS1 CR 
10”^ ctss“^ 
(3) 

MOS2 CR 
10”^ ctss”^ 
(4) 

pnCR 

10“3 ctss“l 
(5) 

HRi 

(6) 

HR2 

0) 

Nr 

(8) 

d(”) 

(9) 

GSC 2.2 

Cat. id. 

(10) 

R 

(11) 

B 

(12) 

V 

(13) 

Nr 

(14) 

d(”) 

(15) 

2MASS 

J 

(16) 

H 

(17) 

Ks 

(18) 

Nr 

(19) 

USNOBl.O 
d(”) Cat. id. 

(20) (21) 

1 

J181731.4-120622.2“ 



311.0±15.3 

0.34±0.05 

-0.67±0.05 


4.4 

S3001223315 


10.6 

9.7 


4.5 

7.3 

7.1 

6.9 


4.4 

0778-0542759 

2 

J181746.1-120542.3 

6.1±1.0 

5.8±1.0 

13.1±1.9 

0.55±0.14 

-0.17±0.16 


2.0 

S300122014588 

17.9 



2 

1.9 

15.0 

13.9 

12.6 


2.3 

0779-0530036 















2.8 

14.7 

13.8 

12.2 




3 

J181750.7-120400.4 

4.4±1.0 

2.6±0.9 

12.1±1.9 

0.51±0.14 

-0.62±0.20 


2.2 

S3001220900 

13.0 

14.4 

13.3 


1.5 

11.4 

11.1 

10.9 


1.9 

0779-0530128 

4 

J181750.7-120506.3 

3.5±0.8 

4.2±0.8* 

11.8±1.8 

0.30±0.14 

-0.47±0.22 


1.7 

S300122015025 

15.7 

18.7 

16.3 


1.4 

13.4 

12.6 

12.4 


2.3 

0779-0530127 

5 

JI81806.0-121433.8'’ 

231.4±4.4* 

229.6±4.1* 

673.3±8.3* 

0.23±0.01* 

-0.75±0.01* 


1.6 

S30012201153 


8.3 

7.6 


1.6 

5.5 

5.3 

5.1 


1.6 

0777-0543906 

6 

J181808.4-120824.9 

1.3±0.4 

1.2±0.4 

5.7±0.9 

0.51 ±0.31 

0.42±0.I5 


2.6 

S300122012953 

17.9 




2.2 

15.6 

15.2 

14.7 

0 



7 

J181808.6-120851.2 

2.6±0.5 

2.1 ±0.4 

3.6±0.8* 

1.00±0.43* 

0.53±0.20* 

0 







2.6 

15.2 

13.8 

13.4 


3.2 

0778-0543325 

8 

J181810.0-121048.5‘= 

12.5±0.8 

13.1±0.9 

36.6±1.9* 

0.38±0.05* 

-0.80±0.05* 


2.1 

S3001220I65 


10.0 

9.3 


2.2 

7.7 

7.5 

lA 


2.1 

0778-0543346 

9 

J1818I0.6-120409.3 

2.1 ±0.5 

2.8±0.6 

7.7±1.2 

0.75±0.16 

-0.16±0.17 


1.8 

S300122015589 

14.3 

16.8 

15.3 


1.5 

11.3 

10.2 

9.4 


1.5 

0779-0530455 

10 

J1818I3.1-121503.1 

4.6±0.7 

4.2±0.7 

3.8±0.9 

0.89±0.16 

-0.82±0.23 


4.4 

S30012209459 

17.5 




3.3 

15.3 

14.4 

14.3 


3.5 

0777-0544085 

11 

J1818I5.3-121148.7 

4.6±0.6 

4.5±0.6 

13.4±1.2 

0.72±0.08 

-0.13±0.09 


1.2 

S300122011032 

18.5 




1.3 

15.2 

14.2 

14.0 

0 



12 

J181815.7-120842.8 

1.2±0.4 

1.5±0.4 

2.6±0.7 

1.00±0.3I 

0.34±0.24 

0 






0 





0 



13 

J1818I6.9-1216I7.4 

3.3±0.6 

1.2±0.5 

3.7±0.9* 

1.00±0.39* 

0.08±0.27* 

0 






0 





0 



14 

J181819.3-121616.9 

5.4±0.8 

3.6±0.6 

13.3±1.7* 

0.54±0.12* 

-0.41±0.14* 

0 






0 





0 



15 

J181820.5-121736.3 

1.7±0.5 

0.6±0.3 

4.2±0.9 

1.00±0.12 

0.16±0.21 

0 






0 






1.1 

0777-0544362 

16 

J181822.6-121503.4 

7.3±0.8 

6.6±0.7 

15.2±1.5* 

0.85±0.17* 

0.45±0.08* 

0 






0 





0 



17 

J181823.2-1207I9.4 

1.1±0.4 

2.2±0.4 

4.1±0.8 

0.38±0.33 

0.42±0.18 

0 






0 






2.5 

0778-0543806 

18 

J181827.3-120613.4 

3.0±0.5 

3.5±0.5* 

11.9±1.2 

0.45±0.09 

-0.67±0.11 


0.9 

S300122014289 

18.1 




0.1 

15.2 

14.1 

13.8 


0.2 

0778-0544019 

19 

J181830.3-121359.1 

4.1 ±0.5 

3.5±0.5 

10.4±1.0 

0.61±0.08 

-0.37±0.11 


1.6 

S30012209935 

17.2 




1.6 

15.2 

14.3 

14.0 


1.6 

0777-0545005 

20 

J181832.0-121740.3 

1.6±0.5 

1.4±0.5* 

4.3±1.0 

0.60±0.20 

-0.60±0.27 


2.4 

S30012208102 

15.2 

18.9 

17.7 

2 

1.4 

13.8 

12.9 

13.3 

2 

2.3 

0777-0545104 















2.8 

13.5 

12.6 

13.1 


2.9 

0777-0545111 

21 

J181832.1-121605.0 

3.5±0.6 

2.4±0.6* 

7.6±1.1 

0.54±0.12 

-0.51±0.16 


3.5 

S300122028708 

18.4 



2 

2.6 

14.7 

13.8 

13.5 


3.3 

0777-0545106 















4.0 

14.2 

13.2 

12.8 




22 

J181837.1-120601.7 

1.6±0.4 

1.3±0.4 

8.1±1.0* 

0.09±0.12* 

-1.00±0.10* 


2.6 

S300122014474 

14.1 

16.2 

15.6 


1.5 

13.5 

12.9 

12.9 


2.8 

0779-0531112 

23 

J181839.5-120935.3 

1.3±0.3 

1.3±3.3 

5.2±0.9* 

0.60±0.16* 

-0.30±0.19* 


3.8 

S30012201041 

13.5 

14.5 

13.8 

2 

1.6 

14.2 

13.0 

12.3 


3.5 

0778-0544645 















4.0 

12.5 

12.2 

12.1 




24 

3181840.9-120623.7'^ 

55.3±1.8 

59.2±1.9 





0.6 

S3001220959 


9.2 

8.6 


0.6 

6.9 

6.7 

6.6 


0.6 

0778-0544725 

25 

J181844.3-120752.3 

12.1±0.9 

13.5±0.9 

45.6±1.9 

-0.37±0.04 

-0.89±0.06 

0 






0 





0 



26 

J181901.2-120050.8 

4.6±0.9 

3.9±0.8 

16.4±2.2 

0.76±0.12 

-0.06±0.14 


1.0 

S300122017698 

15.2 

18.5 

16.6 


1.2 

12.0 

11.0 

10.5 


1.0 

0779-0532345 

27 

J181903.0-120638.4 

4.3±0.8 

3.2±0.6 

12.2±1.5 

-0.26±0.11 

-0.70±0.27 

2 

2.3 

S3001220971 


10.5 

9.8 

2 

2.4 

8.3 

8.1 

8.0 


2.3 

0778-0546021 









3.4 

S300122014341 

12.0 




4.2 

13.8 

12.9 

12.7 




28 

1181904.9-120918.7'= 


2.8±0.6 

7.3±1.3* 

1.00±0.18* 

0.54±0.13* 

0 







2.0 

14.2 

13.1 

12.8 

0 



29 

J181905.2-120306.5 

2.9±0.7 

5.1±0.9 

9.1±1.7 

0.99±0.09 

-0.07±0.18 

0 






2 

0.9 

15.0 

14.0 

13.5 

2 

1.1 

0779-0532614 















3.8 

15.3 

14.3 

14.0 


2.6 

0779-0532625 

30 

J181906.8-1210I8.2 

5.4±0.9 

7.0±0.9 

18.0±1.9* 

0.60±0.19* 

0.60±0.08* 

0 







4.5 

15.4 

14.6 

14.1 


3.0 

0778-0546278 

31 

J181909.9-121753.4 

5.7±1.1 

8.2±1.3 

12.8±2.1 

1.00±0.05 

-0.23±0.17 


1.0 

S300122028430 

17.3 




1.3 

15.1 

14.2 

13.9 


1.1 

0777-0547481 


“ This X-ray source is associated to the 06V((f)) star HD 167834 (BD-12° 4969). We quote the count rates obtained as an extended source detection. The X-ray properties quoted are from 
EPIC-pn data of the April observation only, as this source was not in the field of the other data sets. 

^ HD 167971. The large difference between the count rates quoted for this star in this table and the observed ones given in Table|4|are explained by the exposure map correction applied during 
the source detection procedure. 

BD -12° 4982 (09.51 or III, see text). 

HD 168112 : the EPIC-pn count rate is not quoted as a CCD gap crosses the sources region. The analysis described by De Becker et al. <2004b> was performed with different screening criteria 
than for the current analysis to deal with the EPIC-pn data. 

° We note that even if no optical counterpart was found in the catalogues used for the correlation, a star appears clearly on the DSS image centered on this position. 

* This value is possibly affected by a CCD gap located at <10 arcsec, but is nevertheless quoted. 


(N 
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Wavelength (A) 

Fig. 7. Optical normalized spectra obtained between 4000 and 4900 A respectively for the sources #1 (HD 167834, top spectrum) 
and #8 (BD -12° 4982, bottom spectrum). 


of those objects that have quality flags A or B for the mea¬ 
surements of all three individual near-IR magnitudes. We used 
the March 2003 update of the colour transformations, initially 
derived by Carpenter <200 1> and available on the 2MASS web¬ 
site^, to convert the J — H and H — Kg colours to the ho¬ 
mogenized photometric system introduced by Bessell & Brett 
jl988t . We thus exclude all objects that have colours and mag¬ 
nitudes that are either subject to large uncertainties or are only 
upper limits due to non-detections. Using the extinction law of 
Rieke & Lebofsky <1985> and assuming Ry = Ay /E{B — 
V) = 3.1 (Barbon et al. l2000l . we also show the reddening 
vector for E(B — V) = 1.02 as was found to be appropriate 
for NGC 6604 (Barbon et al. l2000l . 

There are three broad groups that appear in Fig.|3 The 
hrst group consists of five sources that are clearly associated 
with early-type stars: #1 (HD 167834), #5 (HD 167971), #8 
(BD -12° 4982), #24 (HD 168112) and #27 (BD -12° 4994). 
In the second group, we hnd objects around H — K ~ 0.3 
and J — E[ ^ 1.0. These objects have colours that are 
consistent with slightly reddened late type (early G to late 
K) main-sequence or giant stars. Finally, a small group of 
objects (#9, #26 and one of the two counterparts of the three 
sources #21, #23 and #29) have H — K > 0.4. Assuming 
that they are affected by the same reddening as NGC 6608 
(E{B — V) = 1.02), these objects would be associated with 
very late (mostly late M) main sequence stars. If these sources 
were normal main sequence stars belonging to NGC 6604, it 
seems rather unlikely that we would be able to detect their 
X-ray emission. An alternative possibility could be that they 
are X-ray bright pre-main sequence stars with a moderate IR 

http : //WWW.ipac.caltech. edu/2mass/index.html 


excess. In fact, their dereddened infrared colours are in broad 
agreement with the intrinsic colours of classical TTauri stars 
(Meyer et al. ll997> . 


Of course, a fraction of the X-ray selected objects might 
in fact be foreground stars or background sources. NGC 6604 
(III = 18.26°, bii = 1.69°) lies very close to the Galactic 
plane and the total Galactic neutral hydrogen column density 
along this direction must therefore be quite large. This should 
produce a substantial absorption of X-ray photons from extra- 
galactic background sources. Because of the Galactic coordi¬ 
nates of NGC 6604, estimating this column density accurately 
is a very difficult task. For instance, the DIRBE/IRAS extinc¬ 
tion maps provided by Schlegel et al. <1998t yield a good es¬ 
timate of the total Galactic E{B — V) at Galactic latitudes 
above |6ii| > 5°, but are subject to very large uncertainties near 
the Galactic plane. With these limitations in mind, the maps 
of Schlegel et al. (I1998t suggest E{B — U) ~ 3.5 ± 0.5 for 
NGC 6604 (corresponding to Nn ~ 2 x 10^^ cm“^). Assuming 
that extragalactic background sources have a power-law spec¬ 
trum with a photon index of 1.4, and are subject to a total in¬ 
terstellar absorption of 2 x 10^^ cm“^, the detection limits ~ 
4.0 X10“^ and 1.5x 10“^cts s“^ for the pn andMOS detectors 
translate into unabsorbed fluxes of 2.1 x 10“^^ ergcm“^ s“^ 
and 6.0 x 10“^"^ ergcm“^ s“^ in the 0.5 - 2.0keV and 2.0 - 
lOkeV band respectively. From the logW - log S' relation of 
Giacconi et al. (IMD, we would expect to first order to de¬ 
tect around 8-10 extragalactic sources over the EPIC held of 
view. We emphasize that these objects should be detected as 
rather hard sources. Hence, several of the hard sources quoted 
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H -K 


Fig. 8. JHK colour-colour diagram of the 2MASS counter¬ 
parts of the X-ray sources in the EPIC field of view around 
NGC 6604. The heavy solid line yields the intrinsic near- 
IR colours of main sequence stars following Bessell & Brett 
(GUD, whereas the reddening vector is illustrated for Ay = 
3.16. Open circles and triangles stand for EPIC sources having 
respectivley a single or two 2MASS counterpart(s). The dotted 
straight line yields the locus of dereddened colours of classical 
TTauri stars according to Meyer et al. \\991l . 


the EWs of the Hel A 4388 line (0.36 A) and of the Hell A 
4686 line (0.31 A), we derive a supergiant luminosity class, 
whilst the EWs of the Si IV A 4088 line (0.54 A) and of the 
He 1 A 4143 line (0.29 A) point to a giant luminosity class. We 
note however that the spectrum of this star is very similar to 
that of the 09.7Iab star p,Nor (HD 149038) given by Walborn 
& Eitzpatrick ( I1990> . Einally, for the source #14, possibly as¬ 
sociated to Cl* NGC 6604 ED 61 (following the nomenclature 
proposed by Eorbes & DuPuv [T978t even if no optical counter¬ 
part is quoted in Table|^ we derive a G6 spectral type, with an 
uncertainty of about 2 spectral types. We note that this optical 
counterpart is however located at about 6 arcsec of the X-ray 
source #14. Two possible 2MASS infrared counterparts are ap¬ 
proximately located at the position, but the precision of their 
photometry is very poor. 


Source #1 : 2—T MEKAL model 



0.5 1 2 5 

channel energy (keV) 


in Tab le|3 that lack an optical counterpart may actually be as¬ 
sociated with AGN. 

Optical spectra of the counterparts of EPIC sources #1, #8 
and #14 were obtained on June 22, 2004 with the EMMI in¬ 
strument mounted on ESO’s 3.5 m New Technology Telescope 
(NTT) at La Silla. The EMMI instrument was used in the 
RILD low dispersion spectroscopic mode with grism # 5 (600 
groovesmm“^) providing a wavelength coverage from about 
3800 to 7020 A with a spectral resolution of 5.0 A (FWHM of 
the He-Ar lines). The slit width was set to 1 arcsec and the ex¬ 
posure times were 1 min for the two bright sources #1 and #4, 
and 40 min for the counterpart of source #14. The data were 
reduced in the standard way using the long context of the MI¬ 
DAS package. 

The optical spectrum of source #1 (HD 167834) is given in 
the top part of Eig.0 We obtain equivalent widths (EW) of 0.36 
and 0.69 A respectively for the He I A 4471 and He II A 4542 
lines. We estimate that the typical error on the estimate of EWs 
is about 5-10%. According to the classification criterion given 
by Mathys J1988> . and considering that the N III AA 4634-41 
lines are in emission and the He II A 4686 line is in strong ab¬ 
sorption, we derive an 06V((f)) spectral type. We note also that 
the C III A 5696 line is in weak emission. In the case of source 
#8 (BD -12° 4982, bottom spectrum of Eig.Q, the ratio of the 
EW of the He I A 4471 and He II A 4542 lines leads to an 09.5 
spectral type. This spectral type is in excellent agreement with 
that given by Barbon et al. J2000I their source #32). Eollowing 


Fig. 9. EPIC-pn spectrum of source #1 
(HD 167834) obtained in April 2002, fitted with a 
wabsisM*wind* (mekal+mekal ) model. The ISM 
absorption column was frozen at 0.73 x 10^^ cm ^ 


6.3. X-ray properties 


Using the count rates obtained in the three energy bands men¬ 
tioned in Sect. l6.ll we obtained hardness ratios (see columns 6 
and 7 of Table|^ defined respectively as 


HRi 


M-S 
M + S 


HR2 = 


H-M 
H + M' 


We find that some sources present rather extreme hardness ra¬ 
tios, suggesting very hard X-ray spectra. Sources #7, #12, #13, 
#15, #28, #29 and #31 have HRi values larger than 0.90 as 
shown in column 6 of Tablel^for EPIC-pn. We note that high 
hardness ratios are also obtained from EPIC-MOS instruments 
for these sources. Some of these sources may be extragalctic 
sources as discussed above. 

We selected the brightest sources, and performed a spec¬ 
tral analysis of our EPIC data from the April and September 
observations. 
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Table 7. Parameters for models fitted to the spectra of sources #1 between 0.4 and lO.OkeV, and #8 and #25 between 0.4 and 
S.OkeV. The data are from the EPIC-pn observation for source #1, and from the three EPIC instruments for sources #8 and 
#25. The parameters have the same meaning as in Tables|2] and |3 Eor sources #1 and #8, two absorption columns were used 
accounting respectively for the ISM and the local absorption (same ionized wind absorption model as for HD 167971, left as a 
free parameter). The ISM columns were frozen at 0.73 x 10^2 and 0.55 X 10^2 cm 2 respectively for sources #1 and #8. In the 
case of source #25, a unique absorption column (ISM + local) was used and left as a free parameter. The error bars represent the 
90% conhdence interval. The upper and lower parts of the table are respectively devoted to the purely thermal models and to the 
model including a power law. The observed fluxes quoted in the last column are affected by large errors due to the large relative 
error on the normalization parameters (Normi and Norm 2 ), and no signihcant variability between the April and September 
observations can be claimed on the basis of these values. 


Source 

Obs. 

Nh 

(1022 cm"2) 

kTi 

(keV) 

Normi 

kT2 

(keV) 

Norm2 

xl (d.o.f.) 

Obs. Flux 
(ergcm-2 s”*) 

#1 

April 

o.9iJii 

0.19g:?? 

6.20}®6f X 10-2 

1 qo2.02 

^••JOl.07 

1.20j;il X 10"* 

0.97 (78) 

8.0 X 10-^* 

#8 

April 


o.7ig:i 

0.26|:?i X 10“* 

- 

- 

0.73 (66) 

6.8 X 10-1^ 


September 

o.69g:!| 

Q T1 0.80 

' ^0.64 

0.30i:f| X 10“* 

- 

- 

0.85 (178) 

7.9 X 10-1^ 

#25 

April 

(A 'y^o.84 

0.23g:?i 

1.12*:«1 X 10-* 

- 

- 

0.63 (95) 

7.9 X 10-1^ 


September 

n '7/^0.93 
/O0.59 

o.24g:fi 

0.90^i^ X 10-* 

- 

- 

0.87 (133) 

6.1 X 10-1^ 



Nh 

(1022 cm"2) 

kT 

(keV) 

Normi 

r 

Norm2 

xl (d.o.f.) 

Obs. Flux 

(ergcm-2 ^-i^ 

#1 

April 

0.20j:g^ 

0.58g:t| 

1.43o:t7 X 10-* 

2.98l;g 

4.04^2^?® X 10-^ 

1.16(78) 

7.9 X 10-** 


Source #1 is likely associated with the 06V((f)) star 
HD 167834 (BD-12° 4969). This X-ray source was de¬ 
tected by the EINSTEIN satellite with a count rate 
of 0.012 ±0.004cts s“^, leading to a flux* of about 
3.4 X 10~^*ergcm“2 (Grillo et al. ll992t . If we take into 
account the exposure map correction, this source is the second 
brightest X-ray emitter in NGC 6604, after HD 167971 and be¬ 
fore HD 168112. 

Source #8 has a single counterpart in the three catalogues 
used for the source identification, and is most probably asso¬ 
ciated to the 09.51 or III star BD -12° 4982. No optical or in¬ 
frared counterpart was found for source #25. No information 
on the X-ray emission of these two objects has been found in 
the literature. 

In the case of #1, a 50" radius circular region was used 
to extract the spectrum. The EPIC-pn spectrum was binned 
to have a minimum of 9 counts per channel. Eor #8 and #25, 
the spectra were extracted within a circular region with a ra¬ 
dius of 40". Eor the three sources, the background was ex¬ 
tracted within an annulus centered on the source and covering 
the same area as the source region. Gaps and bad columns were 
excluded by rejecting properly adjusted rectangular boxes. The 
spectra were fitted following the same approach as described 
in Sect. 13.1.^ Eor sources #8 and #25 we obtained very sim¬ 
ilar fitting results from the EPIC-MOS and EPIC-pn spectra. 
All spectra from sources #8 and #25 were binned to reach a 
minimum of 5 counts per channel. The main results from the 


The flux was estimated assuming a thermal model with 
kT = 0.5keV and a neutral absorption column obtained by integrat¬ 
ing a model of the hydrogen particle spatial distribution in the Galaxy. 
The procedure is described in Grillo et al. <19921 . 


spectral analysis of these sources are summarized herebelow 
and in Table0 

- #1. Prom the (B — V) = 0.95 color given by Porbes <20001 
for this star, and the intrinsic color typical for an 06.5 star 
(Mihalas & Binnev n'981> . we followed the same procedure 
as described in Sect. B.l.ll to derive E{B — V) = l .26, yield¬ 
ing an ISM A^h of 0.73 x 1022cm“2. We used the same 
ionized wind absorption model as for HD 167971 to ac¬ 
count for local absorption, and we left it as a free parameter. 
The EPIC-pn spectrum is best fitted with a two-temperature 
thermal mekal model (see Fig.|3, with characteristic tem¬ 
peratures of about 0.2 and 1.4keV. The rather high temper¬ 
ature for the hotter component suggests that the intrinsic X- 
ray emission from the shocks in the stellar wind of the OB 
star (see e.g. Feldmeier et al. nw7t is not the only process 
responsible for the observed X-rays. From this model, we 
derive an observed flux of ^ 8.0 x 10“^* ergcm“2 be¬ 
tween 0.4 and lO.OkeV. Assuming that the distance to the 
star is 1.7kpc (Barbon et al. l2000t . we obtain an Lx cor¬ 
rected for the ISM absorption of 4.4 x 10*^ ergs“^. If we 
consider the color excess derived above in this paragraph, 
the V magnitude from Table|^ the distance to NGC 6604 
provided by Barbon et al. <2000t . and the typical bolomet- 
ric correction for an 06V star (Vacca et al. ll996> . we ob¬ 
tain a bolometric luminosity of about 1.7 x 10^®ergs“^. 
This yields an Lx/Lboi ratio of 2.6x 10“®. Considering 
the X-ray luminosity expected from the empirical relation 
of Berghoferet al. <1997t . e.g. ~ 2.8 x 10*2 ergs“^, we ob¬ 
tain an X-ray luminosity excess of about 16! We mention 
that if we use the typical bolometric luminosity for a single 
06V star (i.e. 9.8 x 10*®ergs“^, Howarth & Prinia ll91^ . 
HD 167834 appears a factor 30 overluminous in X-rays. 
































16 


M. De Becker et al.: An XMM-Newton observation of NGC 6604 


Source ffS : 1 —T MEKAL model 


Source #25 : 1-T MEKAL model 



0.5 1 2 5 

channel energy (keV) 



o 

cz 

O 


LD 

LO 

I 

O 

' 0.5 1 2 5 

channel energy (keV) 



Fig. 10. Combined EPIC-MOS1 and EPIC-MOS2 spectra of the September 2002 observation, fitted with a 
wabsisM*wind*mekal model for source #8 (left) and a wabs*mekal model for source #25 (right). In the case of 
source #8, the ISM absorption column was frozen at 0.55 x 10^^ cm 


The spectrum of this source is also compatible with a model 
including a thermal component and a power law, but the 
quality of the fit is somewhat poorer than that obtained 
for the purely thermal model (see Table^}. The possible 
explanations to produce additional X-rays associated to a 
plasma temperature of a few 10^ K could be an interaction 
between the winds of two OB stars in a binary system (see 
e.g. Stevens et al. ll992> . or possibly a magnetic confine¬ 
ment scenario as proposed by Babel & Montmerle J 19971 . 
However, the huge luminosity excess can more probably be 
explained by the binary scenario. Higher quality data are 
needed to discuss the physical origin of this X-ray emis¬ 
sion. 

- #8. We adopted the same approach as for source #1 to es¬ 
tablish the value of the ISM hydrogen column density. With 
a (B — V) = 0.64 (Porbes l2000L we derive an ISM A^h of 
0.55 X 10^^ cm“^. A second absorption column was left as 
a free parameter to account for the absorption by the wind 
material. The best fit is obtained between 0.4 and 5.0keV 
with a single temperature mekal model, with a charac¬ 
teristic temperature of about 0.7 keV. We do not find any 
significant difference in the model parameters obtained for 
our two observations. We obtain an oberved flux between 
0.4 and 5.0keV of 7-8 x lO”^'^ ergcm“^ s“^. The fact that 
this source presents a rather soft spectrum is consistent with 
the rather modest hardness ratios reported in Table|6l The 
combined fit of EPIC-MOS spectra is shown in EigllOl 

- #25. As no optical counterpart was found for this source, 
we were not able to determine a priori an ISM hydrogen 
column density in its direction. So we used a unique absorp¬ 
tion column which was left as a free parameter. The spec¬ 
tral analysis of this source between 0.4 and 5.0keV reveals 
a very soft spectrum fitted with a one-temperature model 
with a kT of about 0.2 keV, which is very consistent with 
the HRi and HR 2 values quoted in Table|6l We did not 
find any indication for a significant variability of the model 
parameters between our April and September observations. 


Eigure cni shows the EPIC-MOS spectra fitted simultane¬ 
ously with the 1-T model. 

6.4. X-ray variability 

We ran the edetect_chain metatask on the three EPIC 
data sets simultaneously, separately for the two observations, 
and we carefully inspected the images in both cases to check 
whether the sources listed in Table|^ are detected at both 
epochs. Most of the apparent variability observed for the X-ray 
sources can be explained by the fact that the effective exposure 
times of our two pointings are significantly different (~ 6 ks 
in April and 11 ks in September for EPIC-pn, and ~ 9 ks 
in April and ~13ks in September for EPIC-MOS), and to 
some extent can also be due to gaps or bad columns. However, 
it appears that the sources #19 and #21 are only detected in 
September 2002, whilst sources #15, #28 and #31 are detected 
only in April 2002. This strong variability suggests a flaring 
behaviour. We extracted light curves to investigate their vari¬ 
ability with more details, but these sources are too faint to pro¬ 
vide conclusive results. Source #31 presents apparently a suf¬ 
ficiently high count rate to give a valuable light curve, but we 
note that this source lies very close to the border of the EPIC 
held. Consequently, the observed count rate is about a factor 
of 3 lower than the exposure corrected value quoted in Table|^ 
Among these targets, only #15, #28 and #31 present a very high 
hardness ratio HRi. A comparison of the count rates obtained 
from the two separate observations reveals no significant vari¬ 
ability for the other X-ray sources detected with XMM-Newton. 

We Anally searched for a short time scale variability in the 
cases of the X-ray sources discussed in Section lOl i.e. #1, #8, 
and #25. We extracted light curves with time bins ranging from 
100 s to 1000 s following the same approach as used in the case 
of HD 167971 (see Sect. l3.2.Tll . Variability tests (chi-square, 
pov-test) performed on these light curves did not reveal any sig¬ 
nificant variability during our two XMM-Newton observations 
for any of the three EPIC instruments. 
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6.5. Star formation in NGC 6604? 

As mentioned in the previous subsections, several indica¬ 
tors (possible flaring activity, position in the JHKs dia¬ 
gram...) suggest some of the X-ray emitters in NGC 6604 
may be candidate PMS stars (Feigelson & Montmerle l 1 999> . 
This fact is at first sight contradictory to the statement of 
Barbon et al. J2000t that pre-ZAMS objects are not detected in 
NGC 6604. However, this lack of detection could be partly ex¬ 
plained by the fact that the study of these authors was limited 
to an area of only 2.1 x 3.3 arcmin^ centered on HD 167971. 
Most of the X-ray sources we detect are indeed located outside 
this area. 

The fact that our investigation is only able to reveal 
the brightest X-ray emitters harboured by NGC 6604 could 
strongly bias our investigation of PMS activity in the cluster. 
Let us consider the faintest X-ray sources of our list. Assuming 
a simple thermal model with a kT of about 0.5 keV®, we ob¬ 
tain an unabsorbed flux between 0.4 and lO.OkeV of about 
1.0 X 10~^‘*ergcm“^ s“^. For a distance to the cluster of about 
1.7 kpc (Barbon et al. l2000t . this unabsorbed flux leads to a lu¬ 
minosity for the faintest source detected in our data of about 
3.5 X 10^° erg s“^. This value is of the order of the highest typ¬ 
ical luminosities expected for PMS stars when no flare is occur- 
ing (Feigelson & Montmerle l 1 999t . This suggests that our data 
only allow a detection of the more luminous PMS candidates 
within NGC 6604, and possibly during flaring stages, i.e. at 
their highest X-ray emission state. Deeper X-ray observations 
are needed to investigate the PMS populations in NGC 6604. 

Unlike the situation in NGC 6383 (Rauw et al. l2003t and 
NGC 6231 (Sana et al. l2005l . the X-ray sources in NGC 6604 
reveal no concentration around the more massive cluster mem¬ 
bers. In NGC 6383 and NGC 6231, the fainter X-ray sources 
are mainly associated to PMS stars, thereby suggesting a close 
relationship between the formation of low-mass and high-mass 
stars. In NGC 6604, the rather sparse spatial distribution of the 
members of the cluster already noted in the visible domain, 
show that this is not a mass-segregated open cluster. Following 
the A^-body simulations performed by Bonnell & Bate J2002> . 
gas accretion should indeed have led to a strong contraction 
of the cluster responsible for a significant mass segregation. A 
possible explanation for the absence of notable mass segrega¬ 
tion in open clusters was proposed by Vine & Bonnell (I2003t . 
According to these authors, the strong stellar winds of O and B 
stars in young open clusters result in a significant gas removal 
from the core, where the OB stars are preferentially formed. 
Their simulations show that the changes in the gravitational 
potential of the core as a result of this gas expulsion is able 
to break the mass segregation {core dissolution), allowing stars 
to drift away from the core. The apparently different dynamical 
stages between NGC 6604 and NGC 6383 is not unexpected re¬ 
garding their ages, i.e. respectively 5 ± 2 (Barbon et al. l2000t 
and 1.7 ± 0.4Myr (FitzGerald et al. llOT^ . 


® With an absorbing interstellar column of about 0.6 x 10^^ cm 
i.e. an approximate mean value of those used for the spectral analysis 
of HD 168112 and HD 167971. 


7. Summary and conclusions 

The analysis of our XMM-Newton data from HD 167971 
reveals that the soft part of the X-ray spectrum is thermal 
with typical plasma temperatures of about 0.2 to 0.8 keV, 
i.e. 2 X 10® to 9 x 10® K. The harder part of the spectrum is 
fitted equally well with a thermal or a power law component. 
On the one hand, if the higher energy emission is thermal, 
the plasma temperatures (~ 2.3 x 10^ to 4.6 x 10^ K) are 
compatible with those expected for a shock heated plasma in 
a wind-wind collision zone, where the winds have reached 
speeds close to their terminal velocities. If we consider the 
short period of the close binary system, such velocities are 
not expected at the position of the collision zone of the winds 
of the 05-8V components. Provided the three components 
are indeed physically connected, this argues in favour of an 
interaction with the third more distant companion. On the other 
hand, if the hard energy emission component is non-thermal, 
the photon index is rather high (~ 3), which is reminiscent of 
the case of other non-thermal radio emitters like 9 Sgr (Rauw 
et al.|2003 or HD 168112 (De Becker et al. l2004bt . The X-ray 
luminosity exceeds the value expected for the ‘canonical’ 
Lx/Lhoi relation by a factor 4. This excess is larger than that 
found for HD 168112 (De Becker et al. l2004bt . Finally, we 
report on a weak decrease of the X-ray flux between our two 
observations separated by about five months possibly due to an 
eclipse of the 05-8V H- 05-8V close binary system. 

By considering several results from recent observations 
of massive stars with X-ray satellites {XMM-Newton) and 
radio telescopes (e.g. VLA), we note that the non-thermal 
radio emitters are possibly not the best candidates to display 
a non-thermal emission in the bandpass of current X-ray 
observatories like XMM-Newton or Chandra. In a scenario 
where non-thermal radio emission requires a wind-wind 
collision in a binary system, wide binaries produce too hard 
a thermal spectrum to allow the detection of any putative 
hard non-thermal contribution. On the contrary, short period 
binaries are not expected to produce an observable synchrotron 
radio spectrum because of the huge opacity of the wind 
material for radio photons. They are however likely to produce 
a softer thermal X-ray spectrum than long period binaries, 
offering the possibility to unveil the expected hard non-thermal 
X-ray emission component. In other words, the simultaneous 
detection of a non-thermal emission in the radio and X-ray 
domains seems rather unlikely. More X-ray and radio observa¬ 
tions of massive binaries (O and WR), whatever the thermal 
or non-thermal nature of their radio emission, are needed to 
interpret the non-thermal emission from massive stars. 

Finally, our investigation of the X-ray emission from 
NGC 6604 led to the detection of 29 sources in addition to 
HD 168112 and HD 167971. The cross-correlation of the posi¬ 
tion of these sources with the GSC, the 2MASS and the USNO 
catalogues revealed that most of them have at least one optical 
or infrared counterpart. We describe a more complete spectral 
analysis in the case of the three brightest objects of our source 
list. Two of them are associated to BD -12° 4982 (09.51 or III) 
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and HD 167834 (06.5V((f))) respectively, and the third one has 
no known optical or infrared counterpart. Some indicators of 
stellar formation activity suggest a few sources of our list may 
be PMS objects. Unfortunately most of the sources detected 
in NGC 6604 are too faint to be investigated in details in the 
context of this study. Finally, the X-ray emitters of NGC 6604 
detected with XMM-Newton are not concentrated close to the 
most massive stars of the cluster. This rather sparse spatial dis¬ 
tribution contrasts with that observed in other open clusters like 
NGC 6383 and NGC 6231 where young stellar objects detected 
in X-rays are mainly concentrated close to the most massive 
stars. 
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